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High resolution spectroscopy of the BCD galaxy Haro 15: 
II. Chemodynamics. 
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ABSTRACT 

We present a detailed study of the physical properties of the nebular material in four 
star-forming knots of the blue compact dwarf galaxy Haro 15. Using long-slit and 
echelle spectroscopy obtained at Las Campanas Observatory, we study the physical 
conditions (electron density and temperatures), ionic and total chemical abundances 
of several atoms, reddening and ionization structure, for the global flux and for the dif- 
ferent kinematical components. The latter was derived by comparing the oxygen and 
sulphur ionic ratios to their corresponding observed emission line ratios (the 77 and 77' 
plots) in different regions of the galaxy. Applying the direct method or empirical re- 
lationships for abundance determination, we perform a comparative analysis between 
these regions. The similarities found in the ionization structure of the different kine- 
matical components implies that the effective temperatures of the ionizing radiation 
fields are very similar in spite of some small differences in the ionization state of the 
different elements. Therefore the different gaseous kinematical components identified 
in each star forming knot are probably ionized by the same star cluster. However, 
the difference in the ionizing structure of the two knots with knot A showing a lower 
effective temperature than knot B, suggests a different evolutionary stage for them 
consistent with the presence of an older and more evolved stellar population in the 
first. 

Key words: ISM: abundances - Hii regions - galaxies: abundances - galaxies: funda- 
mental parameters - galaxies: starburst - galaxies: individual: Haro 15 - 



1 INTRODUCTION 



Giant Extragalactic Hll Regions (GHllRs) are extended, lu- 
minous objects, which are observed on the discs of spirals 
and in irregular galaxies. GHllRs are formed due to the pres- 
ence of a large number of young and massive stars whose 
ultraviolet flux ionizes the surrounding gas. The behaviour 
of emission lines observed in Hll galaxies and Blue Compact 
Dwarfs (BCDs) spectra resembles those of Giant Hll Re- 
gions. We can, therefore, use the same analysis techniques 
developed for the latter to derive the temperatures, densi- 
ties and chemical composition of the interstellar gas in the 
star formation bursts detected in these low metallicity galax- 
ies (Sargent & Searle 1970; French 1980; Terlevich et al. 
1991). BCD galaxies present strong starbursts easily iden- 
tified through their strong and narrow emission lines, low 
metallicity environments and a complex history of star for- 



mation. These facts make them interesting objects to study 
metallicity effects in galaxies (Kunth & Ostlin 2000). 

As mentioned in Firpo et al. (2011, hereafter Paper I) 
Haro 15 has an absolute magnitude Mb = -20.69, a surface 
brightness /i= 18.56 mag arcsec"^ and a colour B-V= 0.33 
(Cairos et al. 2001). At a distance of 86.6 Mpc (de Vau- 
couleurs et al. 1991), Haro 15 meets the criteria for a Lu- 
minous Compact Blue Galaxy (LCBG) (Hoyos et al. 2004), 
although some authors suggest that it may represent the fi- 
nal outcome of a merger between a dwarf elliptical and a gas 
rich dwarf galaxy or HI cloud (Ostlin 1998; Gumming et al. 
2008a). The interactions taking place during the merging 
process would act as the starburst trigger. 

Shi et al. (2005) derived an overall average oxygen abun- 
dance of 12-|-log(0/H) = 8.33 for the Haro 15 galaxy, based 
on an electron temperature determination of 8330 K. It was, 
however, previously known from Ha imaging (Cairos et al. 
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2001) that the galaxy shows a complex morphology, where 
the star forming region can be split among a large num- 
ber of knots scattered throughout the whole galaxy. Lopez- 
Sanchez & Esteban (2009b) obtained spectroscopic informa- 
tion from individual regions and determined densities, tem- 
peratures and chemical abundances for several knots within 
Haro 15. Oxygen abundances were derived for regions A and 
B (names following Cairos et al., 2001) and found to be 
8.37±0.10 and 8.10±0.06, respectively. 

Esteban & Vi'lchez (1992) introduced the term "chemo- 
dynamics" after merging chemical and kinematical infor- 
mation in their analysis of the Wolf-Rayet ring nebula 
NGC 6888. In their research, these authors combined high 
resolution spectra with observations at high spatial resolu- 
tion for different areas of the nebula, which enabled them 
to analyse the different velocity components of the studied 
region. 

James et al. (2009) studied the BCD galaxy Mrk996 us- 
ing high spectral resolution data obtained with VLT VIMOS 
in its IFU mode. They found that this galaxy shows multiple 
components in its emission line profiles, evident as relatively 
broad and narrow components, and were able to deblend 
them in order to analyse their physical properties separately. 
According to the authors this peculiar BCD galaxy has ex- 
tremely dense gas in its nuclear region, where stellar out- 
flows and shock fronts contribute to generate the relatively 
broad feature. However James and collaborators could not 
find evidence of complex emission line profiles for the au- 
roral lines sensitive to temperature. The [Olll] A 4363 and 
[N ll] A 5755 emission lines, present in their spectra, showed 
simple profiles which could be fitted by a broad Gaussian 
function. 

In a previous work we presented a detailed study of 
the internal kinematics of the nebular material in multiple 
knots of the blue compact dwarf galaxy Haro 15 (Paper I) 
from echelle spectroscopy. In that paper we have performed 
a thorough analysis of its emission lines, including multi- 
ple component fits to the profiles of their strong emission 
lines. Our results have shown that the giant Hll regions of 
Haro 15 present a complex structure in all their emission 
profiles, detected both in recombination and forbidden lines. 
The emission lines of the brightest region, Haro 15 A, can 
be split in at least two strong narrow components plus a 
broad one. Although the observed regions tend to follow the 
galaxy kinematics, the components of knot A have relative 
velocities that are too large to be explained by galactic rota- 
tion. Almost all knots follow the relation between luminosity 
and velocity dispersion found for virialized systems, either 
when considering single profile fittings or the strong narrow 
components in more complex fits. The one component fits 
show a relatively fiatter slope. 

Generally, the flux calibration of echelle data is difficult 
to check and the weak auroral emission lines do not have 
enough signal-to-noise ratio (S/N) to allow the accurate 
measurement necessary for a reliable estimate of the tem- 
perature of the emitting gas. Long-slit spectroscopic data of 
intermediate resolution are more sensitive and more reliably 
calibrated to perform temperature estimates and abundance 
analysis. The combination of echelle and long-slit data pro- 
vide an opportunity to interpret the abundance results in the 
light of the internal kinematics of the region while providing 
important checks on the spectrophotometric calibration. 



The development of Integral Field Unit (IFU) instru- 
ments has provided the spatial coverage requiered to study 
extended galactic or extra galactic star-forming regions (see 
e.g. Relafio et al. 2010; Cairos et al. 2010; Monreal-Ibero 
et al. 2010, 2011; Rosales- Ortega et al. 2010, 2011; Perez- 
Gallego et al. 2010; Garci'a-Benito et al. 2010; Sanchez et al. 
2011; Lopez-Sanchez et al. 2011; Perez-Montero et al. 2011). 
However, slit spectroscopy (with medium or high spectral 
resolution) is a valid option for spectrophotometric analy- 
sis. This occurs when the object is very compact, or even 
extended but with few star-forming knots, and when good 
spatial and spectral resolution together with a simultane- 
ous wide spectral coverage are required (see e.g. Gumming 
et al. 2008b; Firpo et al. 2010, 2011; Hiigele et al. 2006, 
2007, 2008; Hagele 2008; Hagele et al. 2009, 2010, 2011; 
Perez-Montero et al. 2009; Lopez-Sanchez & Esteban 2009b, 
2010a,b; Lopez-Sanchez 2010). 

In this paper, we perform a detailed analysis of the 
physical characteristics of the ionized gas in Haro 15, dis- 
tinguishing among individual components and global line 
fluxes. Electron temperatures and densities were estimated 
and used to derive ionic and total abundances for all the 
different species such as: O, S, Ne and Ar. In what follows, 
we describe the details of the observations and the relevant 
data reduction procedures. Section 3 explains the different 
methods used to determine the physical conditions of the 
ionized gas which led to the calculation of their chemical 
abundances, together with their uncertainties. The results 
for each region from each applied technique are discussed in 
Section 4, and the summary and conclusions are presented 
in Section 5. 



2 OBSERVATIONS AND DATA REDUCTION 

We obtained a long-slit moderate resolution spectrum using 
the Wide-Field CCD (WFCCD) camera (September 28 of 
2005) mounted on the 100-inch du Pont Telescope, Las Cam- 
panas Observatory (LCO). The observations correspond to 
knots B and C of the BCD galaxy Haro 15 (see Fig. 1 of Pa- 
per 1). The TEK5 detector was used covering the wavelength 
range 3800 to 9300 A (centred at Ac = 6550 A). The effec- 
tive slit width was 1" giving a spectral resolution of R ~ 900 
(AA = 7.5A at A 6700 A), as measured from the FWHM of 
the He-Ne-Ar comparison lines taken for wavelength cali- 
bration purposes. Observing conditions were good, with an 
average seeing of 1" and photometric sky conditions. 

Bias and dome flat-field frames were taken at the 
beginning of the night. The images were processed and 
analysed with IRAF^ routines following standard proce- 
dures. The procedure includes the removal of cosmic rays, 
bias-substraction, division by a normalized flat-field and 
wavelength-calibration. To substract the sky emission we 
defined two windows at both sides of the 2D spectrum of 
each observed knot. For knot B the sky substraction is very 
good. The final spectrum shows relatively weak residuals of 
the OH sky telluric emission lines which not affect any im- 
portant emission line from the star-forming region. In the 



^ Image Reduction and Analysis Facility, distributed by NOAO, 
operated by AURA, Inc., under agreement with NSF. 
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case of knot C, we have a very good sky substraction for 
the blue part of the spectrum (A < 7000 A). We do not 
measure any emission line in the red part of this spectrum. 
The goodness of the sky substractions were taken into ac- 
count in the estimated errors of the measurements of the 
emission lines. The standard star EG 131 was observed for 
flux calibration purposes and the exposure time for this flux 
standard star was 180 seconds. The spectra were corrected 
for atmospheric extinction and flux calibrated. 

A certain degree of second order contamination is 
present in our long-slit spectrum at wavelengths larger than 
about 6000 A. We are not able to deconvolve the different 
contributions of each order in the red end of the spectral 
range, thus, the measured fluxes of the emission lines in that 
part of the spectrum may have been slightly overestimated, 
but these contributions are not very important for the line 
ratios. 

Five different regions in Haro 15 (see Fig. 1 of Paper 
I) were observed with high resolution spectroscopy obtained 
using an echelle spectrograph attached to the 100-inch du 
Pont Telescope, Las Campanas Observatory (LCO), in July 
19 and 20 of 2006. The spectral range of the observations 
covers from 3400 to 10000 A. This spectral range guarantees 
the simultaneous measurement of the nebular emission lines 
from [O II] AA 3727,3729 to [S iii] AA 9069,9532 A at both ends 
of the spectrum, in the very same region of the galaxy. Ob- 
serving conditions were good with an average seeing of 1" 
and photometric nights. 

A detailed description about the echelle data and the 
reduction procedure can be found in Paper I. The effective 
slit width and length were 1" and 4" respectively, and the 
spectral resolution achieved in our du Pont Echelle data 
was R~25000 (AA=0.25A at A6OOOA), as measured from 
the FWHM of the Th-Ar comparison lines taken for wave- 
length calibration purposes. The spectra were obtained as 
single exposures of 1800 seconds each. Flux calibration was 
performed by observing, with an exposure time of 1200 s, 
the CALSPEC spectrophotometric standard star FeigellO 
(Bohlin et al. 2001). This star flux is tabulated every 2A, 
which is ideal for calibrating high resolution echelle spec- 
tra. In addition, Th-Ar comparison spectra, milky flats (sky 
flats obtained with a diffuser, during the afternoon) and 
bias frames were taken every night. A journal of observa- 
tions is shown in Table 1. The corrected data were reduced 
with IRAF routines following procedures similar to those 
described in Firpo et al. (2005). It must be noted that data 
obtained for the fainter knot F have not sufficient signal to 
noise ratio as to allow flux measurements with the accuracy 
needed to perform an analysis of the physical conditions of 
the gas and have therefore been excluded from this paper. 



3 RESULTS 

3.1 Line intensities and reddening correction 

The one dimensional spectra of knots B and C extracted 
from the WFCCD long-slit data are shown in Fig. 1. The 
spectral plots also include some of the relevant identified 
emission lines. 

For these long-slit data the emission line fluxes were 
measured using the splot task in IRAF following the pro- 



Table 1. Journal of observations for the knots. The first column 
is the observation mode, the second column is the nomenclature 
used in Paper I, the third column is the observed date, the fourth 
column is the exposure time for the knots, the fifth column is the 



air masses. 



Mode 


knots 


Date 


Exp. time 


sec z 


long-slit 


B,C 


2005 Sep 28 


2x1200 -1- 1x900 


1.5 


echelle 


A 


2006 Jul 19 


1800 


1.2 




B 


2006 Jul 19 


1800 


1.1 




C 


2006 Jul 19 


1800 


1.1 




E 


2006 Jul 20 


1800 


1.2 



cedure described in Hagele et al. (2006). We used two meth- 
ods to integrate the emission line flux: (i) In the case of an 
isolated line or two blended and unresolved lines, the in- 
tensity was calculated integrating between two points given 
by the position of the local continuum placed by eye; (ii) 
if two lines are blended, but they can be resolved, we have 
used a multiple Gaussian fit procedure to estimate individ- 
ual fluxes. Both procedures are well detailed in Hagele et 
al. (2006; 2008, hereafter H06 and H08, respectively). In the 
case of the echelle data we have also followed the procedure 
described in (i) to measure the flux of a given line including 
all the different kinematical components. This flux will be 
referred to as "global" in this work. 

The statistical errors associated with the observed emis- 
sion fluxes have been calculated using the expression 



0-; 



a,N^^^[l + EW/{NA)]''/'' 



following the procedure described in Perez-Montero & Diaz 
(2003) where ai is the error in the observed line flux, ac 
represents the standard deviation in a box near the mea- 
sured emission line and stands for the error in the contin- 
uum placement, N is the number of pixels used in the mea- 
surement of the line flux, EW is the line equivalent width 
expressed in A, and A is the wavelength dispersion in A per 
pixel (Gonzalez-Delgado et al. 1994). 

Following the iterative procedure presented in Firpo 
et al. (2010), we measured the fluxes and estimated the er- 
rors associated with each kinematical component identified 
in the high resolution echelle data for the strongest emis- 
sion lines (see Paper I). For those lines with S/N too low to 
perform a self-consistent fitting of the different kinematical 
components present in their emission profiles, the iterative 
process used to deconvolve the different components did not 
yield meaningful results. For these weak emission lines we 
used the fits found for the strongest emission lines in Paper 
I with similar ionization degree as the initial guess as in- 
put for the task ngaussf it of IRAF. With these we fit the 
emission profiles of the weak lines fixing the centroids and 
widths of the corresponding initial approximation, allowing 
only the profiles amplitudes to vary. 

We assumed a two ionization zone scheme to select the 
appropriate initial approximation to fit the weak lines: the 
low ionization zone where [Ol], [On], [Nil], [Sll], [Slll] and 
[Ar III] forbidden lines are originated, and the high ionization 
zone where the helium recombination lines and the [O ill] 
and [Neiii] forbidden lines are emitted. We used the self- 
consistent solution found in Paper I for Ha, [N 11] or [S 11] 



4 G. F. Hdgele et al. 





Si 

0-) CO 


' S 


& 




2.5e-15 


=7 C^l 








„ 2e^l5 


_ 








< 


- 






/ 


'a 1.5e-15 


- 








1 


- 




:£■ 




e le-15 


- 




IOIIIJ4363 
14686 




5e-16 


-U 


yiij 


inl 


4\ 





^ 1 




Wv^vJUju, 




^^'^V'^^^M^A^ 



4000 



5000 



6000 7000 

Wavelength ( A) 



8000 



9000 



3e-16 



2e-16 



S le-16 




Mf^Yrt^^ 




4000 



6000 7000 

Wavelength ( A) 



8000 



9000 



Figure 1. Long-slit spectra of knots B and C of Haro 15 (upper and lower panel, respectively). 



(depending on the minimization of the fitting errors) in the 
low ionization zone, and [O ill] in the high ionization one. 

Upper panels of Figs. 2 and .3 show the [Olll] A 5007 A 
emission lines in the Flux- Velocity plane with their kinemat- 
ical components overlapped for knots A and B, respectively. 
In the lower panels of these figures, we show examples of the 
fittings performed for He I A 5876 A and [O ill] A 4363 A for 
knots A and B, respectively, using as the initial approxima- 
tions those shown in the upper panels. We can appreciate 
that this procedure for the weak emission lines gives very 
good fittings with small residuals. 

The existence of an underlying stellar population is sug- 
gested by the detection of absorption features that depress 
the Balmer emission lines in the long-slit and echelle spec- 
tra. Therefore, the errors introduced by this effect in the 
measurement of line intensities were minimized by defining a 
pseudo-continuum at the base of the hydrogen emission lines 
(see H06). The presence of the wings of the absorption lines 
imply that, even though we have used a pseudo-continuum, 
there is still a fraction of the emitted fiux that we are not 



able to measure accurately (see discussion in Diaz 1988). 
This fraction is not the same for all lines, nor are the ratios 
between the absorbed fractions and the emission. In H06 
it was estimated that the difference in Balmer lines fiuxes 
calculated using the defined pseudo-continuum or a multi- 
Gaussian fit to the absorption and emission components lies 
within the errors. At any rate, for the Balmer emission lines 
we have doubled the estimated error, a;, as a conservative 
approach to consider the uncertainties introduced by the 
presence of the underlying stellar population. 

The reddening coefficient [c(II/3)] was calculated assum- 
ing the galactic extinction law of Miller & Mathews (1972) 
with i?v=3.2. We obtained a value for c(H/?) in each case by 
performing a least-square fit to the observed ratio between 
F(As) and F(H/3), where B denotes the different Balmer 
lines, to their theoretical values computed by Storey & Hum- 
mer (1995) using an iterative method to estimate Ue and Te 
in each case (Hagele 2008). Whenever available, initial Ue 
and Te guesses were taken from measured [S ii] AA 6717,6731 
and [Olll] AA 4363,4959,5007 line fluxes. We have considered 
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Figure 2. [O in] A 5007 A (upper panel) and He I A 5876 A (lower 
panel) emission lines of knot A in the Flux- Velocity plane with 
their kinematical decomposition and the sum (hatched area) of 
the different components superposed. 




Haro 15 B 
[O m] ?.5007 



6400 
Velocity fkm s "1 




6400 

Velocity (km s ) 



Figure 3. [O III] A 5007 A (upper panel) and [Olll] A 4363 A (lower 
panel) emission lines of knot B in the Flux- Velocity plane with 
their kinematical decomposition and the sum (hatched area) of 
the different components superposed. 



He equal to n([Sn]). Owing to the large error introduced by 
the presence of the underlying stellar population, only the 
four strongest Balmer emission lines (Ha, H/3, H7 and US) 
have been taken into account. 

The emission line fluxes of the observed knots of Haro 15 
are listed in Table 2 for the long-slit data, and in Tables 
3-6 for the echelle data. Each table lists the reddening cor- 
rected emission lines ratios for each measurement, either for 
global or kinematical component, together with the redden- 
ing constant and its error taken as the uncertainties of the 
least square flt, and the reddening corrected H/3 intensity. 
First two columns share the same information in Tables 2-6. 
Column 1 lists the wavelength and identification of the mea- 
sured lines. The adopted reddening curve, f(A), normalized 
to H/3, is given in column 2. The following columns show the 
equivalent widths (EW) in A, the reddening corrected emis- 
sion line intensities relative to H/3, and their corresponding 
errors obtained by propagating in quadrature the observa- 
tional errors in the emission line fluxes and the reddening 
constant uncertainties for different knot/component accord- 
ingly. We have not taken into account errors in the theoret- 
ical intensities since they are lower than the observational 
ones. 



As the [S in] A 9532 A line is outside the spectral range 
of the long-slit data, all the physical parameters of knot 
B that depend on this emission line were calculated using 
the theoretical ratio between this line and [S ill] A 9069A, 
I(9532)Ri2.44xI(9069) (Osterbrock 1989). In the case of knot 
C, we were not able to measure the emission line flux of 
[S III] A 9069A with acceptable accuracy. Thus, we do not in- 
clude this line either in the table or in our calculations. For 
this knot, it was impossible to measure the [Nil] A 6548 A 
emission line since it is blended with Ha. We corrected the 
Ha flux for the contamination by this nitrogen line esti- 
mating this contribution from [N 11] A 6584 A using the the 
theoretical ratio, 1(6584) « 2.9 • 1(6548) (Osterbrock 1989). 

We can appreciate in the echelle data that the 
[S 11] A 6731 A emission line is affected by OH sky telluric 
emission lines present in the same spectral range due to the 
galaxy redshift. We used the standard star Feige 110 (sdOB) 
as telluric star since this star has no stellar lines in this 
spectral range. Using IRAF telluric routing, we have cor- 
rected the echelle spectra by this effect. The following step 
was to resort to the ngaussf it task using the fitting solu- 
tion for [Sll] A 6717 A as the initial approximation to fit the 
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Table 2. Relative reddening corrected line intensities [F(H/3)= 
B and C. 



:/(-ff/3)=10000] for the global measure of the long-slit spectrum of knots 









B 






C 




A (A) 


f(A) 


-EW 

(A) 


/(A) 


Error 

{%) 


-EW 

(A) 


/(A) Error 

(%) 


3727 [On]" 


0.271 


109.1 


11940 ± 120 


1.0 


45.3 


25400 + 2100 


S.3 


3835 H9 


0.246 


3.8 


470 ± 20 


5.0 








3868 [Neiii] 


0.238 


34.0 


4490 ± 10 


0.2 








3889 Hei+H8 


0.233 


9.4 


1260 ± 70 


5.1 








3970 [Nelll]+He 


0.215 


21.1 


2535 ± 4 


0.2 








4068 [S II] 


0.195 


0.7 


90 ± 20 


19.6 








4102 H(5 


0.188 


21.1 


2400 ± 30 


1.3 








4340 H7 


0.142 


41.5 


4480 ± 20 


0.4 


6.0 


3850 + 660 17.1 


4363 [O III] 


0.138 


6.3 


690 ± 40 


6.0 








4471 He I 


0.106 


2.8 


310 ± 20 


5.2 


3.8 


2840 + 570 19.9 


4658 [Feiii] 


0.053 


0.3 


36 ± 6 


17.8 








4686 He 11 


0.045 


1.5 


150 ± 20 


12.1 


1.6 


1260 + 330 26.1 


4713 [Ariv]+Hei 


0.038 


0.6 


64 + 7 


11.4 








4740 [Ariv] 


0.031 


0.9 


90 + 10 


14.1 








4861 H/3 


0.000 


114.1 


10000 + 30 


0.3 


14.6 


10000 + 540 


5.4 


4959 [O III] 


-0.024 


221.3 


21220 ± 80 


0.4 


8.3 


6680 + 210 


3.1 


5007 [Oiii] 


-0.035 


539.1 


51340 ± 90 


0.2 


27.9 


22040 + 470 


2.1 


5876 He I 


-0.209 


22.5 


1240 + 30 


2.3 








6300 [Oi] 


-0.276 


6.8 


370 + 20 


4.3 


2.5 


1540 + 300 19.1 


6312 [Siii] 


-0.278 


2.3 


120 + 20 


15.5 








6364 [Oi] 


-0.285 


2.2 


120 + 10 


9.6 








6548 [Nil] 


-0.311 


8.3 


420 + 20 


4.8 








6563 Ha 


-0.313 


565.8 


27620 ± 40 


0.1 


52.6 


28500 + 820 


2.9 


6584 [Nil] 


-0.316 


23.4 


1180 + 20 


1.8 


5.8 


3230 + 370 11.5 


6678 He I 


-0.329 


5.5 


260 + 10 


4.7 








6717 [S 11] 


-0.334 


32.1 


1450 + 20 


1.7 


10.2 


5210 + 610 11.8 


6731 [Sii] 


-0.336 


13.7 


620 + 50 


8.0 


5.1 


2610 + 490 18.9 


7065 He I 


-0.377 


6.3 


250 + 20 


7.2 










7136 [Arm] 


-0.385 


19.4 


710 + 20 


3.0 










7281 Hei' 


-0.402 


1.1 


42 + 5 


11.9 










7319 [Oii]'= 


-0.406 


4.0 


148 + 6 


3.9 










7330 [On]'' 


-0.407 


3.1 


116 + 4 


3.2 










7751 [Arm] 


-0.451 


5.8 


190 + 10 


7.0 










9069 [Siii] 


-0.561 


69.9 


1200 + 120 


9.8 










I(H/3)(ergs-^cm 
c(H/3) 


-'') 




4.49 xlO"^* 
0.29 + 0.01 






0.11 xlO"^* 
0.24 + 0.12 




" [0 11] AA 3726 + 3729; " 
= [Oii] AA 7318 + 7320; <* 


possibly blend with an unknown line 
On] AA 7330 + 7331. 







[S 11] A 6731 A line profile fixing the centroids and widths of 
the components, allowing only to fit the amplitudes. 

The emission lines with relative flux errors higher 
than 45% were not taken into account for the calcula- 
tions. The exception is [Olll] A 4363 A auroral emission line 
in knot B, which, albeit faint, is essential to estimate the 
[O III] line temperature. We were not able to measure the 
[O 11] AA 7319,7330 A auroral emission lines for any echelle 
spectra. In knots C and E, the [S ill] A 9069,9532 A lines were 
not detected due to the low signal-to-noise of these spec- 
tra in this spectral range. In these knots the [Olll] A 4959 A 
emission line fluxes were calculated from the [O ill] A 5007 A 
line using the theoretical ratio, 1(5007) ^ 3x 1(4959) (Oster- 
brock 1989). In knot A, the weak H5 emission line was not 
considered for calculations. For knot E, the narrow 2 com- 
ponent and the global measure present c(H/3) values slightly 
lower than zero but still compatible with this value within 
the observational errors. In this case we considered c(H/3) as 
zero. 



3.2 Physical conditions of the ionized gas 

We have studied the physical conditions of the emitting 
medium in four star-forming knots of Haro 15 using long-slit 
(only for knots B and C) and echelle data. With the echelle 
data, we have analysed the conditions in each knot from the 
global and kinematical decomposition measurements follow- 
ing the results derived in Paper I. In what follows we will 
explain the different methods used to estimate the electron 
temperatures (Te) and densities [Ne~n([Sn])] (see discus- 
sion in Firpo 2011, and references therein). 

(i) When the weak auroral lines needed to implement the di- 
rect method could be measured we have derived the physi- 
cal conditions using the direct method as described in H08 
and Hagele et al. (2011, hereafter Hll). This procedure is 
based on the flve-level statistical equilibrium atom approx- 
imation in the task temden, of the software package IRAF 
(De Robertis et al. 1987; Shaw & Dufour 1995). We have 
taken the uncertainties associated with the measurement of 
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Table 3. Relative reddening corrected line intensities [F{Hl3)=I{Hf})- 
nents of the echelle spectrum of knot A. 



:10000] for the global measure and different kinematical compo- 









global 






narrow 1 






narrow 2 






broad 




A (A) 


f(A) 


-EW 


7(A) 


Error 


-EW 


/(A) 


Error 


-EW 


/(A) 


Error 


-EW 


/(A) 


Error 






(A) 




(%) 


(A) 




(%) 


(A) 




(%) 


(A) 




(%) 


3727 [On]" 


0.271 


12.8 


30070 ± 1700 


5.7 


6.7 


34390 ± 1570 


4.6 


2.5 


37850 ± 4180 


11.0 


1.9 


33500 ± 5980 


17.9 


3868 [Neiii] 


0.238 


1.3 


2900 ± 290 


10.1 




















3888 H8+Hei 


0.233 


0.8 


1280 ± 200 


15.8 




















3970 [Nelll]+He 


0.215 


1.7 


1940 ± 360 


18.7 




















4102 US 


0.188 


5.4 


5340 ± 1300 


24.4 




















4340 H7 


0.142 


6.2 


5140 ± 280 


5.4 


2.2 


4590 ± 190 


4.1 


0.8 


4210 ± 550 


13.0 


0.4 


2800 ± 940 


33.6 


4471 He I 


0.106 


0.4 


470 ± 90 


19.6 




















4861 H^ 


0.000 


12.9 


10000 ± 340 


3.4 


6.2 


10000 ± 160 


1.6 


2.9 


10000 ± 290 


2.9 


2.1 


10000 ± 590 


5.9 


4959 [Oiii] 


-0.024 


6.3 


6300 ± 120 


2.0 


2.3 


4500 ± 170 


3.8 


2.3 


9240 ± 440 


4.8 


1.4 


7800 ± 940 


12.0 


5007 [Oiii] 


-0.035 


22.4 


19400 ± 220 


1.1 


8.3 


14190 ± 160 


1.1 


7.5 


26530 ± 320 


1.2 


5.4 


26720 ± 980 


3.7 


5876 He i 


-0.209 


1.7 


990 ± 60 


6.1 


0.6 


690 ± 80 


10.9 


0.7 


1350 ± 140 


10.5 


0.5 


1540 ± 350 


22.5 


6300 [Oi] 


-0.276 


1.7 


760 ± 40 


5.7 




















6548 [Nil] 


-0.311 


4.2 


1860 ± 70 


3.6 


2.1 


1910 ± 80 


4.3 


1.1 


1440 ± 110 


7.3 


0.7 


1650 ± 450 


27.3 


6563 Ha 


-0.313 


74.0 


30100 ± 220 


0.7 


31.8 


28500 ± 140 


0.5 


22.1 


28500 ± 210 


0.7 


12.9 


28500 ± 660 


2.3 


6584 [Nil] 


-0.316 


13.3 


6010 ± 170 


2.9 


5.9 


5570 ± 130 


2.4 


3.5 


4810 ± 200 


4.2 


3.2 


7430 ± 690 


9.3 


6678 He I 


-0.329 


1.0 


420 ± 50 


11.4 




















6717 [Sii] 


-0.334 


8.2 


3800 ± 120 


3.2 


3.1 


2940 ± 100 


3.5 


1.4 


1780 ± 130 


7.3 


3.8 


8530 ± 840 


9.9 


6731 [Sii] 


-0.336 


8.3 


3250 ± 210 


6.4 


1.9 


1890 ± 120 


6.5 


1.1 


1430 ± 130 


9.2 


2.2 


5250 ± 670 


12.7 


7065 He I 


-0.377 


0.6 


270 ± 90 


32.9 




















7136 [Arm] 


-0.385 


2.1 


860 ± 50 


6.1 


0.7 


580 ± 80 


14.5 


0.5 


520 ± 80 


16.1 


0.7 


1440 ± 380 


26.4 


9069 [Siii] 


-0.561 


4.4 


2170 ± 140 


6.3 


1.6 


1800 ± 210 


11.5 


1.5 


1850 ± 250 


13.6 


0.5 


1260 ± 390 


31.0 


9532 [Siii] 


-0.592 


5.6 


3110 ± 250 


8.1 


1.1 


1450 ± 250 


17.0 


2.2 


2860 ± 330 


11.6 


1.7 


4630 ± 1180 


25.4 


I(H/3)(ergs-^cm-^) 




1.56 X lO"^* 






0.85 X 10-1-* 






0.39x10-1* 






0.28 X 10-1"' 




c(H/3) 






0.22 ± 0.04 






0.11 ± 0.03 






0.68 ± 0.05 






0.37 ± 0.12 





' [0 11] AA 3726 -I- 3729 



the emission-line fluxes and the reddening correction as error 
sources, and we have propagated them through our calcu- 
lations. The adopted coUisional atomic coefficients are the 
same used in H08. 

(ii) When some of the weak auroral lines are not detected in the 
spectrum of an observed region, and the line temperatures 
cannot be computed using the direct method, we have used 
relations between temperatures derived using photoioniza- 
tion models (Perez-Montero & Diaz 2003, 2005), or empir- 
ical and semi-empirical relations (Hagele et al. 2006; Diaz 
et al. 2007; Pilyugin 2007). Perez-Montero & Diaz (2003, 
2005) examined some relationships between temperatures 
used for ionic abundance determinations, adapting these re- 
lations to the latest photoionization model results at that 
moment and using the most recent atomic coefficients. We 
have to stress that the errors of the line temperatures derived 
using these model dependent relations are formal errors cal- 
culated from the measured line intensity errors applying er- 
ror propagation formulae, without assigning any error to the 
temperature calibration itself. We have estimated line tem- 
peratures following the method described in HOG and H08 
which takes advantage of the temperatures derived using 
the direct method in combination with either model based 
or empirical relations. 

(iii) When no auroral emission line was detected in the spec- 
trum, we have used empirical methods to estimate the elec- 
tron temperatures: 

• If the strong emission lines of [Sii], [Sin], [On], and 
[O in] are available in the spectrum of the studied region, 
we can use the relation between the empirical parameter 



SO23 (=823/023; defined by Diaz & Perez-Montero 2000) 
and the [S in] line temperature developed by Diaz et al. 
(2007, hereafter D07): 

ie([5iii]) = 0.596 - 0.283 log 5*023 + 0.199 (log 5023)^ 

This relationship was derived for high metallicity regions. 
However, from Fig. 9 of D07, we can appreciate that this 
relation is also valid for the H n galaxies in the upper end 
of the SO23 distribution for this kind of objects (where 
their dispersion is the lowest) and whose derived values 
of SO23 are in the validity range of the relation. Fortu- 
nately this is the case for the SO23 empirical parameter 
as measured in the star-forming knots of Haro 15. 

We used this Te([S in]) to estimate Te([0 in]) using the 
relation between these two line temperatures derived by 
H06: 

ie([S'iii]) = (1.19 ± 0.08)te([Oiii]) - (0.32 ± 0.10) 

The temperature errors have been calculated from the 
measured line intensity errors, taking also into account 
the temperature calibration errors. 

• If the [S in] emission lines are not detected with 
enough quality, we can use the relation between the [N 11] 
nebular-to-auroral line intensities as a function of the in- 
tensities of the strong nebular oxygen lines for the deter- 
mination of the [Nil] electron temperature (t2) applying 
the empirical relation by Pilyugin (2007, hereafter P07): 

1.111 
*^ ^ logQivn - 0.892 - 0.144 ■ log ti + 0.023 • ti 
where Qjvn is the approximation to the [N n] line intensi- 
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Table 4. Relative reddening corrected line intensities [i<'{H/9)=7(H/3)=10000] for the global measure and different kinematical compo- 
nents of the echcUe spectrum of knot B. 



A (A) 



f(A) 



-EW 

(A) 



global 
7(A) 



Error 

{%) 



-EW 

(A) 



7(A) 



Error -EW 

(%) (A) 



broad 

7(A) Error 

(%) 



3727 [On]" 
3835 H9 
3868 [Neiii] 
3889 He i-|-H8 
3968 [Neiii]-|-H7 
3970 [Nelll]-|-He 
4102 H(5 
4340 H7 
4363 [O III] 
4471 He I 
4713 [Ariv]+Hei 

4861 np 

4959 [Oiii] 
5007 [Oiii] 
5876 He I 
6300 [OI] 
6312 [Siii] 
6548 [Nil] 
6563 Ha 
6584 [Nil] 
6678 He: 
6717 [Sii] 
6731 [Sii] 
7065 He I 
7136 [Arm] 
7751 [Arm] 
9069 [Siii] 

9532 [Siii] 

I(H/3)(ergs^-' cm" 

c(H/3) 

° [0 11] AA 3726 - 



0.271 
0.246 
0.238 
0.233 
0.216 
0.215 
0.188 
0.142 
0.138 
0.106 
0.038 
0.000 
-0.024 
-0.035 
-0.209 
-0.276 
-0.278 
-0.311 
-0.313 
-0.316 
-0.329 
-0.334 
-0.336 
-0.377 
-0.385 
-0.451 
-0.561 
-0.592 



20.4 
1.0 
12.4 
5.1 
5.2 
3.8 
9.7 
37.8 
4.4 
3.0 
2.6 
81.3 
130.0 
470.0 
12.4 
2.7 
1.8 
2.4 
358.0 
7.7 
4.8 
10.6 
8.4 
4.9 
11.0 
3.2 
13.9 
21.3 



9050 ± 240 

360 ± 20 

5590 ± 180 

1570 ± 40 

1500 ± 30 

1060 ± 60 

2770 ± 120 

4830 ± 120 

780 ± 130 

470 ± 20 

340 ± 30 

10000 ± 100 

22710 ± 180 

69810 ± 210 

1053 ± 9 

181 ± 8 

130 ± 30 

185 ± 5 

28620 ± 40 

620 ± 10 

320 ± 3 

730 ± 10 

590 ± 20 

290 ± 6 

710 ± 9 

200 ± 10 

1190 ± 110 

2510 ± 140 

3.16 X 10"^* 

0.13 ± 0.01 



2.6 
4.4 
3.3 
2.6 
2.0 
5.3 
4.3 
2.5 
16.3 
3.6 
8.4 
1.0 
0.8 
0.3 
0.9 
4.4 
20.2 
2.6 
0.1 
1.9 
1.1 
1.3 
3.8 
1.9 
1.3 
6.4 
9.3 
5.6 



8.1 24100 ± 850 3.5 8.0 5880 ± 310 5.3 
1.6 4210 ± 650 15.4 8.4 5850 ± 270 4.6 



0.5 
1.2 
1.4 
5.6 
0.4 



14.1 
16.7 
35.0 
1.7 
2.0 
0.7 
1.6 
73.5 
6.0 
0.5 
7.4 
5.2 

4.0 

3.1 
6.9 



1020 ± 380 
1630 ± 240 
2220 ± 340 
5310 ± 160 
470 ± 210 



10000 

14960 

37510 

570 

500 

200 

420 

29350 

1570 

130 

1920 

1410 



± 160 
± 290 
± 800 
± 40 
± 40 
± 40 
± 20 
± 270 
± 50 
± 30 
± 50 
± 50 



980 ± 50 



1010 
2890 



140 
190 



36.9 
14.9 
15.2 
3.0 
45.8 



1.6 
1.9 
2.1 
6.5 
7.5 
22.2 
5.7 
0.9 
3.4 
20.7 
2.4 
3.8 

5.3 

13.5 
6.4 



2.8 
3.4 
5.7 
16.7 
2.4 



45.6 
91.0 
233.8 
10.2 
0.8 
1.0 
1.6 
184.8 
4.0 
4.6 
3.6 
2.6 

6.9 

9.2 
13.1 



1500 ± 160 
1260 ± 100 
2530 ± 110 
4430 ± 60 
720 ± 90 



10000 ± 70 

25570 ± 150 

79550 ± 340 

1240 ± 20 

70 ± 20 

100 ± 30 

160 ± 20 

28010 ± 130 

400 ± 30 

430 ± 20 

360 ± 30 

270 ± 30 

670 ± 20 

1330 ± 90 
2490 ± 130 



TT 



2.40 X 10" 1" 
0.07 ± 0.01 



10.8 
8.2 
4.5 
1.3 
12.5 



0.7 
0.6 
0.4 
1.6 
30.3 
24.1 
13.0 
0.5 
6.2 
3.8 
8.1 
12.6 

3.4 

6.8 
5.0 



') 



0.74 X 10 



0.36 ± 0.02 



-3729 



Table 5. Relative reddening corrected line intensities [F{Hf3)=I{Hf})- 
nents of the echelle spectrum of knot C. 



:10000] for the global measure and different kinematical compo- 



A (A) f(A) 



-EW 

(A) 



global 
7(A) 



Error 

(%) 



-EW 

(A) 



narrow 
7(A) 



Error 

(%) 



-EW 

(A) 



broad 
7(A) 



Error 

(%) 



4861 H/3 


0.000 


3.2 


10000 ± 2190 


21.9 


1.0 


10000 ± 1460 


14.6 


1.1 


10000 ± 2350 


23.5 


4959 [Oiii] 


-0.024 


5.6 


6880 ± 390 


5.7 


3.4 


9820 ± 510 


5.2 


2.0 


5250 ± 610 


11.7 


5007 [Oiii] 


-0.035 


5.6 


20470 ± 1270 


6.2 


3.4 


29370 ± 1580 


5.4 


2.0 


15580 ± 1860 


11.9 


6563 Ha 


-0.313 


18.1 


28500 ± 1150 


4.0 


6.5 


28500 ± 710 


2.5 


8.0 


28500 ± 800 


2.8 


6584 [N 11] 


-0.316 


1.3 


2330 ± 740 


31.7 


0.9 


3850 ± 940 


24.4 


0.4 


1410 ± 590 


41.7 


6717 [S 11] 


-0.334 


2.6 


5020 ± 1590 


31.7 


1.2 


5400 ± 1360 


25.1 


1.4 


5320 ± 1890 


35.6 


6731 [Sii] 


-0.336 


1.8 


2890 ± 1070 


37.0 


0.7 


3190 ± 1300 


40.6 


1.0 


3890 ± 1610 


41.3 



0.08 X lO-i"* 
0.21 ± 0.40 



0.04 X 10-1" 
0.14 ± 0.25 



0.04 X 10-1'' 
0.38 ± 0.41 



I(H/3)(ergs- 
c(H/3) 



ties ratio (= [Nil] AA6548,84/[N 11] A5755). Qjvn depends 
on the strong [0 11] AA 3727,29 A and [O iii] AA 4959,5007 A 
emission lines. As the calibration of Pilyugin's method has 
been derived using data on high metallicity regions and its 
validity could be questioned in the low metallicity regime, 
we confirmed it in those regions were both methods could 
be applied, finding a good agreement. It is therefore possi- 
ble to estimate the line temperature of [N 11] , which appar- 
ently characterizes the low ionization area of the nebula, 



using the Pilyugin's method. D07 found a systematic dif- 
ference of about 500 K between the [N 11] line temperature 
using the Pilyugin's method and [S in] temperatures de- 
rived using their empirical relation. This difference was 
accounted for when these determinations were present in 
our work. 

(iv) When we cannot measure either any auroral emission line, 
or any intense lines to derive the electron temperatures from 
empirical methods, for example in the two weakest star- 
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Table 6. Relative reddening corrected line intensities [F{Hl3)=I{Hf})=10000] for the global measure and different kinematical compo- 
nents of the echcUe spectrum of knot E. 









global 






narrow 1 






narrow 2 




A (A) 


f(A) 


-EW 


7(A) 


Error 


-EW 


7(A) 


Error 


-EW 


7(A) 


Error 






(A) 




(%) 


(A) 




(%) 


(A) 




(%) 


3727 [On]" 


0.271 


8.2 


32110 ± 9750 


30.3 














4861 H/3 


0.000 


7.5 


10000 ± 1530 


15.3 


1.8 


10000 ± 1270 


12.7 


1.9 


10000 ± 1800 


18.2 


4959 [Oiii] 


-0.024 


2.9 


5750 ± 730 


12.7 


2.2 


3440 ± 350 


10.2 


4.4 


6900 ± 380 


5.6 


5007 [O III] 


-0.035 


6.9 


12160 ± 920 


7.6 


2.2 


10160 ± 1050 


10.3 


4.4 


20670 ± 1210 


5.9 


5876 He I 


-0.209 


1.6 


1730 ± 560 


32.4 














6300 [OI] 


-0.276 


2.1 


2020 ± 560 


27.9 














6548 [Nil] 


-0.311 


1.4 


1420 ± 470 


33.0 














6563 Ha 


-0.313 


32.0 


28500 ± 1090 


3.8 


15.9 


28500 ± 640 


2.2 


11.0 


28500 ± 960 


3.4 


6584 [N II] 


-0.316 


5.8 


5270 ± 1160 


22.0 


1.7 


2810 ± 620 


22.0 


3.7 


9030 ± 2220 


24.6 


6717 [S II] 


-0.334 


6.0 


6210 ± 1440 


23.2 


1.4 


2510 ± 510 


20.4 


4.3 


11230 ± 2880 


25.6 


6731 [Sii] 


-0.336 


5.2 


4450 ± 1400 


31.5 


0.8 


1880 ± 790 


42.1 


1.8 


6640 ± 2630 


39.6 


I(H/3)(ergs- 


-icm-2) 




0.18 X 10-1-* 






0.07 X 10-1"' 






0.07x10-14 




c(H/3) 






0.00 






0.59 ± 0.22 






0.00 





" [O II] AA 3726 + 3729 



forming knots, we considered Te([0 ill]) equal to 10* K, since 
it is a typical value for this kind of objects. Theoretical and 
empirical relations were therefore used to derive the other 
line temperatures from the adopted value for Te([0 ill]). 

3.3 Chemical abundance derivation 

To study the global and kinematical component abundances 
in each knot, we have derived the ionic chemical abun- 
dances of different species. We have used the strongest avail- 
able emission lines detected in the analysed long-slit and 
echelle spectra and the task ionic of the STSDAS package 
in IRAF, based on the five- level statistical equilibrium atom 
approximation, as described in H08. 

The total abundances have been derived taking into 
account, when required, the unseen ionization stages of 
each element, resorting to the most widely accepted 
ionization correction factors (ICFs) for each species, 
X/H = ICF(X+')X+7H+ (see Perez-Montero et al. 2007). 
The procedure is detailed in the following subsections. 

3. 3. 1 Helium 

The helium recombination lines arise mainly from pure re- 
combination, although they could have some contribution 
from coUisional excitation and be affected by self-absorption 
(see Olive & Skillman 2001, 2004, for a complete treat- 
ment of these effects). The electron temperature Te([Olll]) 
is taken as representative of the zone where the He emission 
arises since ratios of recombination lines are only weakly 
sensitive to electron temperature. We have used the equa- 
tions given by Olive and Skillman to derive the He^/H"*" 
value, using the theoretical emissivities scaled to H/3 from 
Benjamin et al. (1999) and the expressions for the coUisional 
correction factors from Kingdon & Ferland (1995). To cal- 
culate the abundance of twice ionized helium we have used 
equation (9) from Kunth & Sargent (1983). A summary of 
the equations used to calculate these ionic abundances is 
given in Appendix B of Garcia-Benito (2009). 

It is possible to calculate the abundances 
of once and twice ionized helium using the 



Hei AA 4471,5876,6678,7065 A and Hen A4686A emis- 
sion lines, respectively. When these lines are available, 
we have used them to derive the corresponding ionic 
abundance. As the observed objects have low densities, 
three of the used helium lines have a small dependence with 
optical depth effects, then we have not made any corrections 
for the fluorescence. We have not corrected either for the 
presence of an underlying stellar population. The total 
abundance of He has been found by adding directly the two 
ionic abundances, He/H= (He'''-|-He^"'")/H^. The results 
obtained for each line and the total He abundances, along 
with their corresponding errors are presented in Table 7 
and 8 for the long-slit and echelle data, respectively. These 
tables also include the adopted value for He"'"/H"'' as the 
average value, weighted by the errors, of the abundances 
derived from each He I emission line. 



3.3.2 Ionic and total chemical abundances from forbidden 
lines 

We have calculated the ionic and total abundances of O, 
S, N, Ne, and Ar using the estimated line temperatures 
as described in H08. The Fe ionic abundance has been 
derived using equations (11) and (12) from Izotov et al. 
(1994). In those cases where there are no measurements 
of the [S ll] A 4068 A and [O ll] A 7325 A auroral emission 
lines, and under the assumption of an homogeneous elec- 
tron temperature in the low-ionization zone, we have taken 
the approximation Te([Sll])f^Te([0 ll]) (Perez-Montero & 
Diaz 2003), where the [On] temperature were derived from 
the [O III] ones using the relations given by photoionization 
models (Perez-Montero & Diaz 2003). The nitrogen auro- 
ral line, [Nil] A 5755 A, is unobservable for all the studied 
star-forming knots. Then, under the same assumption, we 
have estimated this temperature using the approximation 
Te([Nll])«Te([Oll]), except in the cases where we had to 
derive the [N ll] temperatures first using Pilyugin's method 
and estimate the [O ll] ones using the same approximation. 

(i) The oxygen ionic abundance ratios, 0+/H+ and 
0^+/H+, have been derived from the [On] AA 3727,29 A 
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and [Olll] AA4959, 5007 A lines, respectively, using for each 
ion its corresponding temperature. At the temperatures 
derived for the observed objects, most of the oxygen is 
in the form of O^ and O^"""; therefore the approximation 
0/H=(0+ + 0^+)/H+ is valid to represent the total oxy- 
gen abundances. 

(ii) The sulphur ionic abundances, S+/H+ and S^"'"/H+, 
have been derived using Te([Sll]) and Te([Slll]), and 
the fluxes of the [Sll] AA 6717,6731 A and the near-IR 
[S III] AA 9069,9532 A emission lines, respectively. A rela- 
tively important contribution from S'^'*" may be expected for 
sulphur, depending on the nebular excitation. Taking into 
account this unseen ionization state, the total sulphur abun- 
dance is calculated using an ICF for S^ + S^""" . A good ap- 
proximation for this ICF in terms of O^/O is the formula by 
Barker (1980), is based on Stasihska (1978) photoionization 
models, with a = 2.5, which gives the best fit to the scarce 
observational data on S''"'' abundances (Perez-Montero et al. 
2006). H08 further reduced the propagated error by writing 
this ICF in terms of the ratio 0^+/0 instead of 0+/0. 

(iii) The nitrogen ionic abundance, N+/H^, has been 
derived from the intensities of the [Nil] AA 6548,6584 
lines. The N/O abundance ratio was calculated as- 
suming that N/0 = N^/0^ and N/H are estimated as 
log(N/H) = log(N/0) + log(0/H). 

(iv) The neon ionic abundance, Ne^"*", has been derived us- 
ing the [Ne ill] A 3868 A emission line, which is the only vis- 
ible and uncontaminated emission line of this atom in the 
spectra. For this ion we have considered the electron tem- 
perature of [Olll], as representative of the high- ionization 
zone (Te([Nelll])f»Te([Olll]); Peimbert & Costero 1969). 
Izotov et al. (2004) point out that this assumption can lead 
to an overestimate of Ne/H in objects with low excitation, 
where the charge transfer between O'^"'' and H'' becomes im- 
portant. Thus, we have used the expression of this ICF given 
by Perez-Montero et al. (2007). Given the high excitation of 
the observed objects, there is no significant differences be- 
tween the derived neon abundances using this ICF and those 
estimated with the classical approximation (II08, Hll). 

(v) The argon ionic abundance, Ar^"*", has been calculated 
from the measured [Arm] A 7136 A emission line assuming 
that Te([Arlll]) «Te([Slll]) (Garnett 1992). In some cases, 
it is possible to measure the [Ariv] A A 471 1,4740 A emis- 
sion lines. However, the first one is usually merged with the 
He I A 4713 A line, and we are not able to separate them. 
Therefore, we use the second and more intense line to cal- 
culate the Ar^"*" abundance. The ionic abundance of Ar'^+ 
has been calculated assuming that Te([Arlv]) ~Te([Olll]). 
The total abundance of argon are hence calculated using 
the ICF(Ar2+) and the ICF(Ar2+ + Ar3+) derived from pho- 
toionization models by Perez-Montero et al. (2007). 

(vi) Finally, for iron we have used the [Fe ill] A 4658 A 
emission line using the electron temperature of [O ill] 
[Te([Felll])«Te([Olll])]. We have taken the ICF(Fe2+) 
from Rodriguez & Rubin (2004). We were only able to mea- 
sure the [Fe III] A 4658 A emission line flux in the long-slit 
spectrum of knot B. 

The ionic abundances with respect to ionized hydrogen 
of the elements heavier than helium, ICFs, total abundances 
and their corresponding errors are also given in Tables 7 and 
8 for long-slit and echelle data, respectively. 



Table 7. Physical properties, and ionic and total chemical abun- 
dances derived from forbidden lines and the helium recombination 
lines for the global measurements of the long-slit data: knots B 
and C. 





B 


C 


n([Sll]) 


100: 


100: 


te([OlIl]) 


1.26±0.03 


1.01±0.02'= 


te([SlIl]) 


1.19±0.12 


0.88±0.13'' = 


te([Oll]) 


1.20±0.04 


1.21±0.01" 


t=([Sii]) 


0.79±0.11 


1.21±0.01" 


te([Nll]) 


1.20*0.04" 


0.93±0.01'= 


12+log(0+/H+) 


7.34±0.05 


7.82±0.06" 


12+log(0^+/H+) 


7.95±0.03 


7.84±0.03= 


12-|-log(0/H) 


8.04±0.03 


8.13±0.04'' 


12+log(S+/H+) 


5.95±0.20 


6.16±0.07" 


12+log.(S^+/H+) 


6.17±0.23 




ICF(S+ -1- S^+) 


1.41±0.03 




12-|-log(S/H) 


6.52±0.22 




log(S/0) 


-1.52±0.22 




12+log(]V+/H+) 


6.18*0.04" 


6.70±0.07'= 


12-|-log(N/H) 


6.88±0.20" 


7.02±0.23'= 


log(N/0) 


-1.16±0.07" 


-1.11±0.09= 


12+log{Ne^+ /H+) 


7.33±0.03 






ICP(iVe^+) 


1.08±0.01 






12-|-log(Ne/H) 


7.36±0.03 






log{Ne/0) 


-0.68±0.04 






12+log(Ar^+ / H+) 


5.65±0.10 






12+los(Ar^+/H+) 


4.62±0.08 






ICF(Ar2++Ar3+) 


1.03±0.01 






12-|-Iog(Ar/H) 


5.70±0.10 






log(Ar/0) 


-2.34±0.11 






12 + log(Fe^+/ir+) 


5.04±0.10 






ICF(Fe^+) 


5.21±0.62 






12-|-log(Fe/H) 


5.76±0.11 






He+/H+{X4i71) 


0.064±0.003 






ire+/ir+(A5876) 


0.096±0.002 






£re+/ir+(A6678) 


0.072±0.003 






£re+/ir+(A7065) 


0.100±0.008 






J?e+/ir + adopted 


0.086±0.015 






He^+/H+(\4686) 


0.0013±0.0002 






(He/H) 


0.087±0.015 







Note. All the values with a superscript are derived us- 
ing models or empirical/semi-empirical relations, and 
the derived errors are underestimated. Densities in 
cm~^ and temperatures in 10* K 

" derived using temperatures predicted by photoion- 
ization models, see text, ''derived using temperatures 
estimated from photoionization models and/or empir- 
ical methods, '^derived using H06 empirical method, 
■^derived using D07 empirical method, "^derived using 
P07 empirical method. 



Table 8. Physical properties, and ionic and total chemical abundances derived from forbidden lines and the helium recombination lines for the global and the different kinematic 
measurements in the echelle data: knots A, B, C and E. 









A 






B 






C 






E 








global 


narrow 1 


narrow 2 


broad 


global 


narrow 


broad 


global 


narrow 


broad 


global 


narrow 1 


narrow 2 




n([Sll]) 


300: 


100: 


200: 


100: 


150±50 


100: 


150: 


400: 


400: 


1400: 


500: 


900: 


400: 




te([OlIl]) 


1.00±0.02'= 


1.08±0.03'= 


1.16±0.04'= 


0.95±0.05'= 


1.20±0.06 


1.22±0.18 


1.12±0.04 


1.0^ 


1.0^ 


1.0^ 


0.97±0.03'= 


1.0^ 


1.0^ 




te([SlIl]) 


0.87±0.05'* 


0.96±0.06'' 


1.06±0.10'* 


0.8±0.12'' 


1.31±0.17 


1.71±0.30 


1.16±0.16 


0.87±0.13'' 


0.87±0.13'' 


0.87±0.13" 


0.83±0.13'*'= 


0.87±0.13'' 


0.87±0.13'' 




te([Oll])" 


0.93±0.01 


1.28±0.01 


1.06±0.01 


1.08±0.03 


1.11±0.02 


1.28±0.10 


1.07±0.01 


0.87±0.01 


0.87±0.01 


0.76±0.01 


0.84±0.01 


0.80±0.01 


0.88±0.01 




t=([Sn])» 


0.93±0.01 


1.28±0.01 


1.06±0.01 


1.08±0.03 


1.11±0.02 


1.28±0.10 


1.07±0.01 


0.87±0.01 


0.87±0.01 


0.76±0.01 


0.84±0.01 


0.80±0.01 


0.88±0.01 




te([Nll]) 


O.QSiO.Ol" 


1.28±0.01° 


i.oeio.oi" 


1.08±0.03" 


1.11±0.02" 


1.2S±0.10'' 


1.07±0.01" 


0.87±0.01'' 


0.87±0.01'' 


0.76±0.01'' 


0.88±0.01' 


0.80±0.01'' 


0.88±0.01'' 




12+log(0+/W+)'' 


7.92±0.06 


7.87±0.05 


7.87±0.10 


8.00±0.15 


7.21±0.08 


7.61±0.21 


7.09±0.08 








8.01±0.19 








12+log(0=+/H-+) 


7.80±0.03'= 


7.55±0.04'= 


7.75±0.05'= 


S.OliO.ll^ 


8.11±0.07 


7.S4±0.19 


8.26±0.05 


7.83±0.03'' 


7.99±0.02'' 


7.71±0.05'' 


7.70±0.08° 


7.53±0.04^ 


7.84±0.02^ 




12+log(0/H)'' 


8.17±0.05 


8.04±0.05 


8.11±0.08 


8.31±0.13 


8.16±0.07 


8.04±0.20 


8.29±0.05 








8.18±0.16'' 








12+log(S+/H+)" 


6.14±0.04 


5.90±0.04 


5.71±0.07 


6.44±0.09 


5.30±0.05 


5.68±0.12 


5.02±0.08 


6.20±0.17 


6.24±0.16 


6.37±0.28 


6.36±0.17 


5.99±0.20 


6.56±0.16 




12+log(S=+/H+) 


6.42±0.18'' 


6.11±0.18'* 


6.20±0.16'' 


6.54±0.2S'' 


5.93±0.12 


5.79±0.13 


6.03±0.13 
















ICF(S+ + 3^+)'' 


1.05±0.01 


1.02±0.01 


l.OSiO.Ol 


1.08±0.02 


1.72±0.02 


1.16±0.01 


2.13±0.07 
















12+Iog(S/H)'' 


6.62±0.13 


6.33±0.13 


6.34±0.14 


6.83±0.21 


6.26±0.11 


6.11±0.13 


6.40±0.12 
















log(S/0)'' 


-1.54±0.14 


-1.71±0.14 


-1.77±0.16 


-1.48±0.25 


-1.90±0.13 


-1.93±0.23 


-1.88±0.13 
















12+log(]V+/H-+) 


6.97±0.03'' 


6.89±0.03'' 


6.78±0.05'' 


7.06±0.10'' 


5.87±0.05'' 


6.25±0.13" 


5.75±0.06'' 


6.57±0.15'' 


6.78±0.12'' 


6.35±0.19'' 


6.92±0.12'= 


6.65±0.11'' 


7.15±0.12'' 




12+log(N/H) 


7.21±0.23'' 


7.06±0.22'' 


7.06±0.22'' 


7.02±0.39° 


6.82±0.32'' 


6.68:" 


6.95±0.28" 








7.09: = 








log(N/0) 


-0.96±0.06'' 


-0.98±0.06" 


-1.09±0.11'' 


-0.94±0.18" 


-1.34±0.09'' 


-I.36i0.25'' 


-1.34±0.10'' 








-1.09±0.23= 








12+loe(Ne^+ /H+ 


) 7.50±0.07'= 








7.49±0.09 


7.34±0.28 


7.62±0.07 
















ICF(Ne^+) 


1.21±0.01'= 








1.07±0.01 


l.lliO.Ol 


1.07±0.01 
















12 + log(Ne/H) 


7.59±0.07'= 








7.52±0.09 


7.39±0.28 


7.65±0.07 
















log(Ne/0) 


-0.58±0.09= 








-0.64±0.11 


-0.65±0.35 


-0.64±0.09 
















12+log{Ar^+ / H + 


6.05±0.19'' 


5.76±0.20'* 


5.63*0.20"^ 


6.36±0.31'' 


5.56±0.11 


5.50±0.14 


5.65±0.14 
















ICF(Ar^+) 


l.lSiO.Ol'' 


1.19±0.01'' 


l.lSiO.Ol'' 


1.13±0.01'' 


1.38±0.02 


l.lliO.Ol 


1.79±0.08 














^ 


12 + log(Ar/H) 


6.11±0.19'' 


5.83±0.20'' 


5.69±0.20'' 


6.41±0.31'' 


5.70±0.11 


5.55±0.14 


5.90±0.14 














log(Ar/0) 


-2.06±0.20'' 


-2.21±0.21'' 


-2.43±0.21'' 


-I.90i0.34'' 


-2.46±0.13 


-2.49±0.24 


-2.39±0.15 














2 


He+ /H+(\4471) 


0.094±0.018° 








0.096±0.003 


















^^ 


He+ /H+(\5876) 


0.072±0.004° 


O.OSliO.OOS" 


O.lOliO.OlO" 


0.111±0.024'' 


0.079±0.001 


0.044±0.001 


0.092±0.001 








0.124±0.040'' 






Ol 


He+/H+(X6678) 


0.109±0.012° 








0.086±0.001 


0.035±0.006 


0.113±0.003 
















He+ /H+{X7065) 


0.107±0.036'' 








0.109±0.005 


















^ 


He+ /H+ adopted 


0.077±0.023° 


0.051±0.005" 


O.lOliO.OlO" 


0.111±0.024" 


0.080±0.017 


0.043±0.006 


0.094±0.014 








0.124±0.040'' 






5r 


Note. All the values with a 


superscript are derived using models or 


empirical/semi-empirica 


relations. 


and the derived errors 


are underestimated. Densities in cm^^ and 


to 


temperatures in 


lO^K 


























£. 


"derived using temperatures 


predicted by photoionization models, see text, ''derived using temperatures 


estimated from photoionization models and/or 


empirical 


methods. 


"^ 


'^derived using H06 empirical 


method, ''derived using D07 empirical method, "^derived using P07 empirical method, -^assumed temperature = lO"* K 






1 

re 

to 
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4 DISCUSSION 

In the following subsections we will discuss the results ob- 
tained from the two different observation modes. 



4.1 Long-slit data 

Ham 15 B 

Four electron temperatures, Te ( [O ill] ) , Te ( [O ll] ) , Te ( [S ill] ) 
and Te([Sll]), were estimated from direct measurements for 
knot B. The good quality of the data translates into small 
rms errors (2%, 3%, 10% and 14%, respectively, for each 
temperature). In the absence of the [Nil] A5755A emission 
line we assumed the relationships between temperatures to 
beTe([Nll]) ^Te([Oll]). 

Knot B shows an electron density value of Ue = 100 
cm~^, typical for this kind of objects. However, the error 
is too large, thus it should be regarded as an estimation of 
the order magnitude of the density. 

In this knot, we have derived the ionic abundances of 
0+, 0^+, S+, S^+, N+, Ne2+, Ar2+, Ar^+ and ¥e^+ using the 
direct method as described in H08. We have also calculated 
the total abundances of O, S, N, Ne, Ar and Fe, and the 
logarithmic ratios of N/O, S/0, Ne/O and Ar/0. 

The oxygen abundance is 8.04±0.03, which is 0.22 times 
the solar value (12+log(0/H) © =8.69; AUende Prieto et al. 
2001). This oxygen abundance value is in very good agree- 
ment, within the observational errors, to that found by 
Lopez-Sanchez & Esteban (2009b), 8.10±0.06. The sulphur 
abundance is 0.15 the solar value (12+log(S/II) = 7.33; 
Grevesse & Sauval 1998), and the nitrogen abundance is 0.09 
of the solar one (12-|-log(N/H) © = 7.93; Holweger 2001). 

The logarithmic N/O ratio found for knot B is 0.28 dex 
lower than the solar value (log(N/0) = -0.88; Asplund 
et al. 2005), and 0.44 dex higher than the typical value shown 
by this kind of objects which present a plateau at log(N/0) 
about -1.6 dex (Pagel et al. 1979a; AUoin et al. 1979) in 
the low metallicity end of the distribution (see discussion 
in Amorin et al. 2010; Perez-Montero et al. 2011, and ref- 
erences therein). The derived value is on the upper range 
of the distribution of the log(N/0) ratio for this kind of ob- 
jects (see Fig. 4). The log(S/0) ratio is -1.52±0.22, in agree- 
ment, within the errors, with the solar value (log(S/0) =- 
1.36; Grevesse & Sauval 1998) (see Fig. 4). The logarithmic 
Ne/O and Ar/O ratios are also consistent with the solar 
ones (log(Ne/0) © =-0.61 dex and log(Ar/0) © =-2.29 dex; 
Grevesse & Sauval 1998) within the errors (see Tables 7 and 
8). 

We have used the intensities of the well detected 
He I AA 4471,5876,6678 and 7065A, and He 11 A 4686 A lines 
to estimated the ionic abundances of He^ and He^""" , respec- 
tively. Then we have derived the helium total abundance us- 
ing these ionic abundances. The adopted value for He^/H"'" 
is 0.086±0.015 and 0.0013±0.0002 for He^+/H+. The He 
total abundance estimated in this region is 0.087±0.015, 
within the typical values found for Hll galaxies (see e.g. H08, 
and references therein). 



Harold C 

In the case of knot C it was impossible to obtain any electron 
temperature using the direct method. In this case, we esti- 
mated the temperatures using the models and the empirical 
relationships, as explained in the previous section. 

The auroral lines are not detected in the knot C 
spectrum due to the low S/N of these data, so we have 
estimated a [N 11] temperature of 9300 ± 100 K following 
the Pilyugin's method (P07). Using the approximation 
Te([Nll])R:i500K-f Te([Slll]) derived by D07, we have es- 
timated Te([Slll]), and derived Te([Olll]) from the em- 
pirical relation fitted by H06. The electron temperature 
Te([Oii]) has been calculated from Te([Oiii]) using the 
relationship based on photoionization models derived by 
Perez-Montero & Diaz (2003), and from the approximation 
Te([Sll]) ~Te([Oll]) we obtained the electron temperature 
of [Sll]. As in the case of knot B, due to the large errors 
involved, the derived electron density of 100 particles per 
cm"^ has to be regarded as a rough estimation. 

In this knot, we were not able to measure the 
[Fe III] A 4658 emission line and therefore no iron abundance 
could be derived. The [Nil] A 6548 A emission line could not 
be measured either and was taken to be 0.33 times that of 
[N III A 6584 A. 



Comparison among knots 

A comparison between the results for both knots from long- 
slit data might suggest, at face value, that the oxygen abun- 
dance could be slightly higher in knot C. However, the 1 a 
errors derived for these oxygen abundances do not take into 
account the errors in the line temperatures introduced by 
uncertainties in the calibrations of the relations based on 
photoionization models. These unaccounted effect will most 
probably make the observed difference negligible. The esti- 
mated N/O ratio for knot C is only marginally larger than 
that for knot B, but still the same within the observational 
errors. 

The temperatures representative of the high-ionization 
zone are higher for knot B than for knot C. The tempera- 
tures for the low-ionization zone show a different behaviour: 
in the case of Te([Oll]) appear to be similar in both knots 
while the Te([Sll]) in knot B is 4200 K lower than in knot 
C. We have to note that the [S 11] temperature was derived 
using the direct method for knot B, but using photoioniza- 
tion models for knot C. This leads to a difference in tem- 
perature among [On] and [Sll] for knot B, although models 
state that we should expect similar values. This behaviour 
was previously found for several objects (see Fig. 2.17 of 
Hagele 2008). Moreover, Perez-Montero et al. (2010) pre- 
sented a self-consistent study in a sample of 10 Hll galaxies 
computing tailored models with the photoionization code 
CLOUDY (Ferland et al. 1998). They found that the elec- 
tron temperature of [S 11] was overestimated by the mod- 
els, with the corresponding underestimate of its abundance, 
pointing to the possible presence of outer shells of diffuse 
gas in these objects that were not taken into account in 
their models. In contrast, the [Nil] temperature shows for 
knot B a higher value than that derived for knot C by about 
2700 K, but again they were estimated using two different 
methods. Only for comparison purposes we have also esti- 
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Figure 4. N/O and S/O as a function of 12-|-log(0/H) (left and right panel, respectively) for knot A and B (open and filled diamonds, 
violet and yellow, respectively), and the objects from H06, H08, and Hll (open circles, turquoise). Open squares are the Hll galaxies 
from Table 13 of H08. The observed knots are marked in the plot with their names plus a letter which denote the observation mode 
("sd": single dispersion) or the kinematic component measure ("g": global, "n": narrow, and "b": broad). The solar values are shown 
with the usual sun symbol: oxygen from AUende Prieto et al. (2001), nitrogen from Holweger (2001) and sulphur from Grevesse & Sauval 
(1998). These values are linked by a solid line with the solar ratios from Asplund ct al. (2005). 



mated this temperature for knot B using Pilyugin's method 
finding Te([Nn]) = 10500 ± lOOK, which is lower than the 
estimated using the models by 1500 K. Comparing those 
[N ll] temperatures derived using Pilyugin's method, knot 
B still shows higher values than knot C by about 1200 K. 

The excitation degree, given by the logarithmic 
0^'''/0^ ratio, is higher in knot B (0.61 dex) than in knot C 
(0.02 dex), in agreement with the estimated values by Lopez- 
Sanchez & Esteban (2009a). The estimated extinction in 
both knots is small, typical of this kind of objects (Hagele 
2008, and references therein). 



4.2 Echelle data 

We analysed the parameters derived from the global and the 
different kinematical component measurements for the four 
knots using the echelle data. 



Haro 15 A 

A large number of recombination and forbidden lines were 
detected in our echelle spectrum, although we were not able 
to measure any auroral emission line. We have therefore used 
the relation between the SO23 parameter and the [Sin] tem- 
perature derived by D07 to estimate Te([Slll]). From this 
temperature we have calculated T^ [O ill] using the relation- 
ship found by H06. Following the classical analysis (see H08), 
we estimated Te ( [0 11] ) , Te ( [S 11] ) and Te ( [N 11] ) using the re- 
lations given by photoionization models. 

The estimated density errors are large, hence these val- 
ues are only considered as an order of magnitude. All the 



electron density determinations are well below the critical 
value for coUisional de-excitation. 

The global measure present a total oxygen abundance 
of 8. 17 ±0.05, 0.3 times the solar value, while for the differ- 
ent kinematical components we obtained 8.04 ±0.05 (nar- 
row 1), 8.11 ±0.08 (narrow 2), and 8.31 ±0.13 (broad), 0.22, 
0.26 and 0.42 times the solar abundance, respectively. There 
seems to be a difference between broad and narrow compo- 
nents above 1 a level, but since the errors in the calibration 
of the relations based on photoionization models are not 
quantitatively established, we can consider the total abun- 
dances in agreement for all components. The derived total 
oxygen abundance for the different kinematical components 
and for the global measure in this knot present the char- 
acteristic low values, within the errors, that are found for 
Hll galaxies: 12±log(0/H) between 7.94 and 8.19 (Terlevich 
et al. 1991; Hoyos & Diaz 2006). 

Taking into account the errors, the total sulphur abun- 
dance of the broad component is about 0.15 dex higher than 
those shown by the narrow ones, but in agreement consid- 
ering errors at 2 a level. On the other hand, the values for 
total nitrogen abundances are in very good agreement for 
all the components. 

The ionic abundances of sulphur, S^/H"*" and S^^/H"*", 
the ionic nitrogen abundance, N^/H^, and the ionic argon 
abundance, Ar^^/H"'", obtained for the broad component 
are higher than the ionic abundances found for the narrow 
components (between 0.17 to 0.61 dex respect to the nar- 
row 1 component, and between 0.28 to 0.73 dex respect to 
the narrow 2 one). Only for the global measure it was possi- 
ble to obtain the ionic and total abundances of [Ne ill] from 
the measure of the neon 3868 A emission line. 
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We have also derived the sulphur and nitrogen abun- 
dances relative to oxygen abundance, S/0 and N/0. If we 
consider the errors, all the S/O derived values are very sim- 
ilar among them (see Fig. 4). The N/O value estimated for 
the narrow 2 component is the lowest, but its difference with 
the other components is negligible. This behaviour of the 
N/O could be an evidence of a common or very similar chem- 
ical evolution for the different kinematical components (see 
discussion in Perez-Montero et al. 2011). It could be also 
indicative that the different kinematical components corre- 
spond to different phases of the same gas, which is in agree- 
ment with the scenarios proposed, for example, by Tenorio- 
Tagle et al. (1996) or Westmoquette et al. (2009). In knot 
A we also obtained an excess of the N/O ratio (see Fig. 4) 
respect to the typical values found for an Hll galaxy. 

For the global measure we can determine the ionic he- 
lium abundance, from the 4471, 5876, 6678, 7065 A emission 
lines. For the three kinematical components we were only 
able to deconvolve the 5876 A emission line. The narrow 1 
component has a lower ionic abundance than the typical val- 
ues found for other similar galaxies, while the values found 
for the other two components and the global measure are in 
the range presented by Hll galaxies (see H08). 

HarolS B 

In this knot we were able to estimate Te([Olll]), and 
Te([Slll]) using the direct method through the deconvo- 
lution of the emission profile for the different kinemati- 
cal components of [O ill] A 4363 A and [S ill] A 6312 A auro- 
ral emission lines. The precision obtained for Te([Olll]) and 
Te([Slll]) are of the order of 5% and 13% for the global mea- 
sure, 15% and 18% for the narrow component, and 4% and 
14% for the broad component, respectively for each temper- 
ature. These are the only two temperatures directly derived, 
even for the global measure. The [O ii] temperature has been 
obtained from Te([Olll]) using the relations given by pho- 
toionization models (Perez-Montero & Diaz 2003). The [Sll] 
and [N ll] temperatures have been estimated using the ap- 
proximation Te ( [S II] ) « Te ( [N II] ) Ri Te ( [O II] ) . 

The errors in the electron densities are large enough as 
to consider our estimate as an order of magnitude. Once 
again, in all cases the electron densities are well below the 
critical value for collisional de-excitation. 

The ionic abundances in the low-ionization zone 
(0+/H+, S+/H+, N+/H+) derived for the narrow compo- 
nent are higher than those shown by the broad one, while the 
ions in the medium- and high-ionization zones (S'^'*'/H+ and 
Ar2+/H+, and 0^+/H+ and Ne2+/H+, respectively) show 
the opposite behaviour, lower values for the narrow compo- 
nent. 

The total oxygen abundance derived for the global mea- 
sure and the kinematical components are in good agree- 
ment taking into account the observational errors, where 
the global measure is approximately equal to the average 
of the derived values for the different kinematical compo- 
nents weighted by luminosity. The values of the total sulphur 
abundances derived for both kinematical components are 
very similar, with the broad component presenting slightly 
larger values than the narrow one. In the case of the to- 
tal nitrogen abundance we derived very similar values for 
the different components, although we have to note that the 



calculated error for the narrow component is too large, and 
therefore this is only a rough estimate of the nitrogen abun- 
dance. The total abundances of Ne derived for the different 
components are in agreement within the errors, and the Ar 
total abundance from the broad component is slightly larger 
(by about 0.07 dex) than that derived for the narrow one. 
As we said above, the temperatures estimated using the re- 
lationships based on photoionization models only take into 
account the errors of the line intensity measurements, but 
without assigning any errors to the relation parameters. We 
could therefore state that all the components present total 
abundances in overall good agreement. 

The S, N, Ne, and Ar abundance relative to oxygen for 
the narrow and broad components, present mainly the same 
values, even comparing these values with those estimated 
for the global measure (see Fig. 4). As in the case of knot 
A, this could be indicative that both components have a 
common or very similar chemical evolution and/or that both 
components represent two different phases of the same gas. 

From the global measure, we have derived the ionic 
abundance of once ionized helium using the same emission 
lines as for knot A. The adopted abundance value is in good 
agreement to what we find in the long-slit spectrum for the 
same knot, and in the range found for other Hll galaxies 
(see II08). We have been able to deconvolve the profiles of 
He I AA 5876 and 6678 A emission lines finding very low val- 
ues for the narrow component, while the broad component 
presents values typical for this class of objects. 

Haro 15 C and Haro 15 E 

In the cases of the star- forming knots C and E, we were 
not able to detect any auroral emission line. Owing to the 
low S/N of these spectra, other some important and rela- 
tively strong emission lines, such as [On] AA3727A, do not 
have enough intensity to be measured, even for the global 
measure. For all the components of knot C and the kine- 
matical components of knot E, we were not able to mea- 
sure the emission lines needed to apply the relationships 
given by D07 and P07 to derive the [S ill] and the [N ll] line 
temperatures, respectively. We assumed that Te([Oiii]) is 
equal to 10*K (as was discussed in §3.2), then we estimated 
the [O ll] temperature in each case using its relation with 
Te([Olll]) given by photoinization models (Perez-Montero 
& Diaz 2003). Te([Sll]) and Te([Nll]) were assumed equal 
to the [On] temperature, and Te([Slll]) was determined us- 
ing its relationship with the temperature of [O ill] found by 
H06. For the global measure of knot E, it was possible to 
measure the [O ll] AA 3727,29 A emission lines. Using these 
emission lines together with the strong [O ill] lines, we ap- 
plied the Pilyugin's method to estimate the electron tem- 
perature of [Nil]. Using the systematic difference of 500 K 
between the temperatures of [N ll] and [S in] found by D07, 
we estimated Te([Slll]). From the [Slll] temperature we de- 
rived the Te([Oiii]) using the relation by II06, then we have 
used the same relations than for the kinematical components 
of knot E to derive Te([Oll]) and Te([Sll]). 

The estimated electron densities for these two star- 
forming knots are rather higher than those found for knots A 
and B. Due to the very large uncertainties present in all the 
cases, these density values are only estimates of the order 
of magnitude of the electron densities. Only for the global 



Haro 15: Physical properties 15 



measure of knot E, it was possible to obtain an estimate of 
the ionic abundance of O^ and, thus, the total abundances 
of oxygen and nitrogen, and the N/0 ratio. 

In contrast with the behaviour found for the 0^^/H"'' 
ionic abundance for the star- forming knots A and B, the 
ionic abundance of the narrow component for knot C is 
higher (~ 0.3 dex) than the one derived for the broad compo- 
nent. This effect could be due to a different ionization degree 
of the kinematical components. In the case of knot E, the 
total oxygen abundance for the global measure is similar to 
the abundances derived for knots A and B, and the helium 
abundance is slightly larger than that found for the other 
knots. 



Comparison among knots 

The total oxygen abundance derived for the global measure 
of knot A is slightly lower than the value obtained by Lopez- 
Sanchez & Esteban (2009b, 8.37 ±0.10). It could be due to 
the different hypotheses considered to determine the elec- 
tron temperatures and the chemical abundances. The total 
oxygen abundance estimated from the global measure for 
knot B is in very good agreement with the value obtained by 
Lopez-Sanchez & Esteban (8.10 ±0.06) for the same knot. 
We have estimated its electron temperature of [O ill] using 
the direct method, and it is in very good agreement with the 
value derived by Lopez-Sanchez & Esteban (12900 ± 700 K). 
Our estimate of Te([On]) is in agreement, although slightly 
lower, with the value derived by Lopez-Sanchez & Esteban 
(12000 ± 500 K) using the relation given by Garnett (1992). 

Recent studies in this galaxy suggest that the ionizing 
star cluster in knot B would have the youngest population 
of the galaxy (~ 100 Myr). This knot presents a blue color 
and high UV emission, indicating a recent star-formation 
activity, supported by the presence of Wolf-Rayet features 
in its spectrum (Lopez-Sanchez & Esteban 2010a). This is in 
agreement with our measurements of the equivalent width 
of the H/3 emission line [EW(H^)] (see Tables 2-6), which 
suggest that knot B is the youngest one since it has the 
highest EW(H/3) (Terlevich et al. 2004). 

On the other hand, we can assume from the analysis 
performed that the total oxygen abundances derived for 
knots A and B are very similar. The ionic oxygen abun- 
dances of O"*"/!!^ derived for knot B using the long-slit spec- 
trum and the global measure of the echelle data are in good 
agreement taking into account the errors involved. Under 
the same considerations, we can assume that the differences 
between the 0^^/H+ ionic abundances and the total oxygen 
abundances derived for the long-slit and for the global mea- 
sure of the echelle spectrum of knot B are negligible. The 
total oxygen abundances derived from both instruments and 
from different components for knots A, B, C, and E are in 
good agreement taking into account the observational er- 
rors. 

The differences found between these star-forming knots 
of Haro 15 are similar (or even smaller) to what is found 
in other works which spatially resolve individual knots that 
belong to Hll galaxies or Blue Compact Dwarf (BCD) galax- 
ies (see e.g. Kehrig et al. 2008; Cairos et al. 2009a; Perez- 
Montero et al. 2009, 2011; Garcia-Benito et al. 2010; Hagele 
et al. 2011). However, in general, these differences were at- 
tributed to the observational uncertainties (pointing errors. 



seeing variations, etc.) or errors associated to the reddening 
correction and flux calibration. The oxygen abundance vari- 
ations were therefore not assumed as statistically significant, 
concluding that there is a possible common chemical evolu- 
tion scenario in all of them. There are even greater differ- 
ences when comparing the estimated abundances of the indi- 
vidual knots with those derived from the integrated spectra 
of the galaxies. For instance, Cairos et al. (2009a) found for 
the integrated spectrum of Mrk 1418 a lower value of direct 
oxygen abundance by about 0.35 dex (equivalent to a factor 
of 2.2) than for knots 1 and 2 of that galaxy. They pointed 
out that while this variation could refiect a real abundance 
difference in different scales (kpc-sized aperture for the in- 
tegrated spectrum and sizes of the order of 100 pc for indi- 
vidual Hll regions), it may also be due to relatively large 
measurement uncertainties for the weak [O ill] auroral emis- 
sion line. Fortunately, as in Hll (where we used a double 
beam long-slit spectrograph), our data are not affected by 
pointing errors and seeing variations, and the other obser- 
vational uncertainties have a second order effect, since the 
echelle spectrograph simultaneously acquire all the observed 
spectral range. Likewise, the errors associated with the mea- 
surements of the weak auroral emission lines are relatively 
small, specially for [O ill] . 

The estimated extinction for all the star-forming knots 
analysed using different instruments and kinematical com- 
ponents (including the global measure) are very similar and 
consistent with low extinction values, found in this kind of 
objects (Hagele 2008, and references therein). 

4.3 Ionization structure 

Vilchez & Pagel (1988) showed that the quotient of 0+/0^+ 
and S"'"/S^''', called "softness parameter" and denoted by 77, 
is intrinsically related to the shape of the ionizing continuum 
and depends only slightly on the geometry. We can perform 
an analysis of the ionization structure only for knots A and B 
since we need information about two consecutive ionization 
stages of two different atoms, reliably detected only in these 
two star-forming knots. 

The purely observational counterpart of the 77 parame- 
ter is the 77' parameter defined by Vilchez & Pagel (1988) as 
the ratio between [Oll]/[Olll] and [Sll]/[Slll]. The relation 
between both parameters is: 

log ri' = log r? - 0.14 te - 0.16 

where te is the electron temperature in units of lO^'K. 

The left panel of Fig. 5 shows the relation between 
log(0+/0^+) and log(S+/S^+) derived for knots A and B. 
The global measures and the different kinematical compo- 
nents of knot A and B are represented with violet and yellow 
diamonds, respectively. In the right panel we show the re- 
lation between log([0 ll]/[Olll]) vs. log ([S 11]/ [S ill]), which 
does not require the explicit knowledge of the involved line 
temperatures for the derivation of ionic ratios and is inde- 
pendent of the methods used to estimate the line tempera- 
tures. The objects studied by H06, H08, and Hll are rep- 
resented with turquoise circles, and the Hll galaxies from 
the literature (see description and references in H08) with 
squares. The Hll galaxies are located in the region corre- 
sponding to high values of the effective temperatures of the 
radiation field of the ionizing star cluster, where logrj are 
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Figure 5. Left panel: log(0+/0^+) vs. log(S+/S'^+) for knots A and B (filled violet and yellow diamonds, respectively), the objects 
studied in H06, H08 and Hll (turquoise circles), and Hll galaxies from the literature as described in H08 (open squares). Diagonals in 
this diagram correspond to constant values oil). Right panel: log([Oll]/[Olll]) vs. log( [S ll] / [S ill] ) , symbols as in the left panel. Diagonals 
in this diagram correspond to constant values of r;'. As in Fig. 4, the letter added to the name of each knot denote the instrument (for 
long-slit spectrum) or the components (for echelle data). 



between -0.35 and 0.2 (see H06). In Fig. 5 diagonal lines rep- 
resent constant values of log r; (left panel) and log 77' (right 
panel) . 



The ionization structures derived from both diagrams, 
rj and rj\ for all the measurements (both instruments and 
different components) of knot B are almost equal, showing 
almost the same values within the observational errors. This 
implies that the effective temperatures of the radiation fields 
that ionize the gas are very similar for different kinematical 
components, then the ionizing star clusters that excite the 
gas that produce the different kinematical components could 
be the same. For knot A we have also derived very similar re- 
sults from both diagrams. In this case the broad kinematical 
component shows an ionization structure slightly displaced 
in the 7;' diagram toward higher effective temperature of the 
ionizing radiation field than those shown by both narrow 
kinematical components (and the global measure) . However, 
all the components are almost in the same zone in the r) di- 
agram if we take into account the large observational errors. 
This fact, points again towards a common star cluster as the 
ionizing source for all the components of knot A. The dif- 
ference in the ionizing structure of these two knots suggests 
a different evolutionary stage, since all the components of 
knot A seem to be located in a region with lower effective 
temperature of the ionizing radiation field than that where 
the components of knot B are placed. This is in agreement 
with the presence of an older and more evolved ionizing star 
cluster in knot A than in knot B, as suggested by Lopez- 
Sanchez & Esteban (2010a). 



4.4 Chemical abundances from empirical 
calibrators 

Different strong line empirical metallicity calibrators are 
commonly used to estimate the oxygen abundances in ob- 
jects for which direct derivation of electron temperatures 
is not feasible. These empirical methods are based on the 
cooling properties of ionized nebulae which ultimately trans- 
late into a relationship between emission line intensities and 
oxygen abundance. These relationships have been widely 
studied in the literature using different strong-line empir- 
ical methods which are based on calibrations of the relative 
intensity of some bright emission lines against the abun- 
dance of relevant ions present in the nebula (see e.g. Garcia- 
Lorenzo et al. 2008; Cairos et al. 2009b; Garcia-Benito 2009; 
Garcia-Benito et al. 2010). 

In Fig. 6, we show the total oxygen abundances as de- 
rived from several strong-line empirical methods, with their 
corresponding errors estimated taking into account the er- 
rors of the line intensities and also the errors given by the 
calibrations of the empirical parameters. We have also plot- 
ted the total oxygen abundances calculated from the elec- 
tron temperatures measured using the direct method, only 
for the long-slit data of knot B, and those estimated from the 
photoionization models and empirical temperature relations. 
The colours of these non-continuous lines (see the electron 
edition of the journal) correspond to each measure or com- 
ponent, thus the black solid and double-dashed-dotted lines 
represent the total oxygen abundances calculated from the 
long-slit data, and turquoise dashed-dotted (green dashed) 
and magenta dotted lines represent the narrow [narrow 1 
(turquoise) and narrow 2 (green) for knot A and E] and 
broad components from the echelle data, respectively. For 
simplicity, the abundance ranges corresponding to the global 
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measure for knot A and B are not shown. This value is plot- 
ted only for knot E since it is the only previous oxygen 
abundance estimation for this knot. 

Among the available strong-line empirical parameters 
to obtain the total oxygen abundance we studied the S23 
(Vilchez & Esteban 1996) using the calibration obtained 
by Perez- Montero & Diaz (2005); O23 (Pagel et al. 1979b), 
and according to the values measured in our data, we used 
the McGaugh (1991) calibration; SO23 ( = 823/023; Diaz & 
Perez-Montero 2000) applying the calibration from Perez- 
Montero & Diaz (2005); O3N2 (AUoin et al. 1979) with 
the calibration established by Pettini & Pagel (2004); S3O3 
(Stasihska 2006); N2 (Storchi-Bergmann et al. 1994) using 
the empirical calibration from Denicolo et al. (2002); and 
Ar303 (Stasihska 2006). 

In the case of knot A, there is a good agreement be- 
tween the total oxygen abundance derived from the global 
measure using the different parameters and the abundance 
derived using the temperatures obtained from models. Simi- 
lar results are found when using empirical temperature rela- 
tions from the global measure. For the narrow 1 kinematical 
component there is a good agreement in the abundance de- 
termination using the empirical parameters: S23, O23, and 
SO23, and that estimated using the temperature derived by 
models and/or empirical temperature relations, while the 
derived abundance using the other empirical parameters 
show small overestimations except the N2 parameter which 
presents the largest difference. The derived abundance for 
the narrow 2 component using the empirical parameters is 
in agreement with that derived using the temperatures ob- 
tained from models and/or empirical relations, except in the 
case of the N2 parameter. The derived abundances from the 
broad component using different empirical parameters are 
very similar to that derived using temperatures obtained 
from models or empirical relations. The exception is the S23 
parameter that yields an abundance slightly higher. 

For the echelle data of knot B, the oxygen abundance 
calculated for the global measure using the empirical pa- 
rameters: O23, O3N2, and N2 shows a similar value to 
that derived using the temperatures obtained from models 
and/or empirical relations. The estimated abundance from 
the other empirical parameters present lower values than the 
abundance derived using the line temperatures. The oxygen 
abundance for the narrow component derived using the em- 
pirical parameters present similar values than that estimated 
using the line temperatures, and for the broad component 
the empirical parameters always present lower values than 
the one determined using the derived temperatures, except 
for the O23 that shows a similar estimation. For the single 
dispersion data, the total oxygen abundance derived from 
empirical parameters are generally in good agreement whit 
that calculated using the temperature obtained using the 
direct method. Only the abundance determined using the 
SO23 empirical parameter is slightly lower than the previ- 
ously estimated. 

For knot C, the total oxygen abundance derived using 
the empirical parameters are in good agreement with the 
oxygen abundance value determined using the temperatures 
derived by models and/or empirical temperatures relations 
for the long-slit data. For knot E, there is a good agreement 
between the abundances derived from the echelle measure- 
ments using the empirical parameters and that estimated 



using the adopted line temperatures, except for the N2 pa- 
rameter using the narrow 2 component. 

We must stress that the different empirical metallicity 
relationships were derived from global measures of complete 
regions, and an analysis or calibration for kinematical com- 
ponents has never been done before. To disentangle the em- 
pirical relation for these two components, it is mandatory to 
observe star-forming regions with high spectral resolution, 
similar to the data presented in our work. Furthermore, even 
better S/N will be needed, in particular to measure the au- 
roral O^ emission lines, to be able to obtain direct mea- 
surements of the total oxygen abundances for the different 
kinematical components. It must be noted that empirical re- 
lations for the estimation of abundances should be applied 
only in a statistical sense and their application to individual 
objects or regions is always questionable. 



5 SUMMARY AND CONCLUSIONS 

We have performed an analysis of the characteristics of 
the ionized gas belonging to the four brightest star-forming 
knots of the BCD galaxy Haro 15. For two of these regions 
we present new low-resolution long-slit data covering from 
3800 to 9300 A obtained with the Wide-Field CCD camera 
mounted on the duPont Telescope at LCO. For all the re- 
gions we have carried out a chemodynamics study using high 
resolution (R ~ 25000) data acquired with the echelle spec- 
trograph mounted at the same telescope. These high resolu- 
tion data were previously presented in a kinematical study 
in Paper I. Using the direct method, and model based and 
empirical temperature relations described in detail in H08, 
we have estimated the electron densities and temperatures, 
ionic and total chemical abundances for different atomic 
species. This analysis was done on the global measurements 
performed to the emission lines of the spectra obtained with 
both instruments, and on the different kinematical compo- 
nents deconvolved in the emission line profiles in the high 
resolution data. 

For all observed knots, and for all the global measure- 
ments and kinematical components, the electron densities 
were found to be well below the critical density for coUisional 
de-excitation, as is well known to occur in the star-forming 
processes belonging to Hll galaxies. 

From the single dispersion spectra, in knot B we 
measured four line temperatures: Te([Olll]), Te([Oll]), 
Te([Slll]) and Te([Sll]) reaching high precision, with rms 
fractional errors of the order of 2%, 3%, 10% and 14%, re- 
spectively. For knot C, from model based and empirical tem- 
perature relations, we estimated the electron temperatures: 
T,([Olll]), Te([On]), T,([Slll]), T,([Slll) and T,([Nll]), al- 
though in these cases the quoted errors for the line tem- 
peratures derived using model based relationships are only 
formal errors calculated from the measured line intensity er- 
rors applying error propagation formula, without assigning 
any error to the temperature calibration itself. 

From the echelle spectra of knot B, we measured 
Te([Olll]) and Te([Slll]) applying the direct method. We 
have reached high precision from the global measure and the 
broad component for the first temperature, with rms frac- 
tional errors of the order of 5 and 4 per cent, respectively, 
and slightly worse, 15 per cent for the narrow component. 
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Figure 6. The total oxygen abundances and their uncertainties for both instruments, and the global measure and different kinematical 
components in each observed knot of Haro 15, as derived using different empirical calibrators. Prom top to bottom panels: knots A, B, 
C, and E. From left to right; the empirical parameters S23, 023> SO23, O3N2, S3O3, N2, and Ar303. Solid black error bars represent 
the long-slit estimate, blue stars: global echelle measure, turquoise up triangles: narrow kinematical component (for knot A and E, filled 
up triangles: narrow 1 component, and filled green squares: narrow 2 component), and magenta down triangles: broad component. The 
solid lines represent the valid range taking into account the 1 a errors of the total oxygen abundances previously derived using the 
direct method, and the non-continuous lines the valid range of the total oxygen abundances derived using photoionization models and/or 
empirical temperature relations (see text). 



We have also obtained good precision for the estimates of 
the [S III] temperature, 13, 18, and 14 per cent from the 
global measure, the narrow and the broad components, re- 
spectively. This is the first time that physical conditions are 
directly estimated for kinematical components of Hll galax- 



Using the estimated values for the electron densities 
and temperatures and a careful and realistic treatment of 
observational errors, we have estimated ionic and total abun- 
dances of O, S, N, Ne, Ar, Fe and He. For the echelle data, 
we have also been able to carry out a chemodynamics anal- 
ysis applying the direct method to the kinematical compo- 
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nent decomposition of the emission line profiles. This kind 
of analysis had never been done until now for this type of 
objects. The obtained total abundances of O, S, N, Ne, and 
Ar are in the typical range found for Hll galaxies. The total 
oxygen abundance derived for knot B is in very good agree- 
ment with the value estimated by Lopez-Sanchez & Esteban 
(2009b). The total oxygen abundances derived for all the 
measurements of the echelle data of this knot are in good 
agreement. For knot A, the derived total oxygen abundance 
from the global measure is lower (~0.2dex) than the value 
found by Lopez-Sanchez & Esteban (2009b). However, tak- 
ing into account all the possible sources of uncertainties we 
find this difference negligible within the errors. Total oxygen 
abundances derived for knots A, B, C, and E are very similar 
among themselves considering the observational errors. The 
N/0 ratios derived from the different components of knot 
A and the global measure of knot E are in very good agree- 
ment, and show an excess with respect to the typical values 
found for Hll galaxies. Knot B seems to also show an excess 
in the N/O ratio, but it is not as evident as in the case of 
knot A. For knots A and B, the relative abundances of N/O, 
S/O, Ne/O, and Ar/0 for the different kinematical compo- 
nents are very similar, these could be evidence for a common 
or very similar chemical evolution for the different kinemati- 
cal components of each knot. It could be also indicative that 
the different kinematical components are different phases of 
the same gas. 

The ionization structure of knots A and B mapped 
through the use of the rj and rj' diagrams shows very similar 
values within the errors for the different components of each 
of these regions. Only the broad component of knot A seems 
to have slightly higher effective temperature (lower value of 
the parameters) than the other components of this knots, 
showing similar values to those of knot B. These similarities 
between the ionization structure of the different kinematical 
components implies that the effective temperatures of the 
ionizing radiation fields are very similar for all the different 
kinematical components, in spite of some small differences 
in the ionization state of different elements. The ionizing 
star clusters that excite the gas belonging to each star- 
forming knot that produce the different kinematical com- 
ponents could therefore be the same. The difference in the 
ionizing structure of these two knots suggests a different evo- 
lutionary stage, with knot A located in a region with lower 
effective temperature than that where knot B are placed. 
This is in agreement with the presence of an older and more 
evolved stellar population in knot A than in knot B, as sug- 
gested by Lopez-Sanchez & Esteban (2010a). 

We have also derived the total oxygen abundances us- 
ing several strong-line empirical parameters. In general, 
the estimated abundances are consistent with the derived 
abundances using the temperatures calculated by the direct 
method and/or by model based and empirical temperature 
relationships. We have to emphasize that all the strong-line 
empirical parameters were derived from measurements of the 
total intensity of the emission lines from entire regions, and 
it has never been done an analysis or calibration for kinemat- 
ical components. Observations with high spectral resolution 
and much better S/N, in particular to measure the auro- 
ral O"'' emission lines, are needed to develop the empirical 
relations for the different kinematical components. 
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